While calculating RKKY interaction in graphene no theorist should do a cut-off 

without cause 
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In our previous work (E. Kogan, Phys. Rev. B 84, 115119 (2011)) we presented calculation 
of RKKY interaction between two magnetic impurities in graphene based on Matsubara Green's 
functions (MGF) in the coordinate - imaginary time representation. Now we present the calculation 
based on MGF in the coordinate - frequency representation. We claim that both approaches have 
an important advantage over those based on zero temperature Green's functions (ZTGF), which are 
very briefly reviewed in the beginning of the present work. The MGF approaches, in distinction to 
the ZTGF approaches, operate only with the convergent integrals from the start to the end of the 
calculation. The coordinate - frequency representation for the MGF turns out to be as convenient 
as the coordinate - imaginary time representation and allows to easily consider the cases of doped 
and gapped graphene. 

PACS numbers: 75.30.Hx;75.10.Lp 
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INTRODUCTION 

RKKY interaction between two magnetic impurities in 
graphene was theoretically studied quite intensely dur- 
ing last several years [ll4l8j, (A terse but precise review 
of the issue one can find in the book by M. Katsnel- 
son 0.) One may ask, why the problem, which is in 
principle so simple (when being treated in the lowest or- 
der of perturbation theory, like it was done in all the 
papers referenced above), was the subject of so many 
publication, using different approaches? The answer to 
this question, as presented below, is connected with the 
fact that a simply written integral is not necessarily a 
simply calculated integral, and in the frame work of all 
zero temperature Green's functions (ZTGF) the integrals 
defining the RKKY interaction in graphene turned out 
to be divergent. Different ZTGF approaches can be thus 
viewed as different ways to obtain finite results from the 
divergent integrals. 

In our previous publication dealing with the subject 
we were using the approach based on Matsubara Green's 
functions (MGF) in the coordinate - imaginary time rep- 
resentation [ll[ . In this work we consider RKKY interac- 
tion in graphene in the framework of approach based on 
MGF in the coordinate - frequency representation, and 
find that this approach is no less convenient than the one 
we used previously. Before we present the approach, we 
give a brief review of the existing ZTGF approaches, in- 
dicating explicitly where the divergent integrals appear. 

If we ignore the spin-orbit coupling, the effective ex- 
change RKKY interaction between the two magnetic im- 
purities with the spins Si and S2, sitting on top of carbon 
atoms at the sites i and j, is 

ij 2 X (i,j)S r S 2 , (1) 
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where J is the contact exchange interaction between 
each of spins and the graphene electrons, and \R is the 



free electrons charge susceptibility. Thus different ap- 
proached to calculation of the RKKY interaction are ac- 
tually different approaches to calculation the susceptibil- 
ity. 

Not to distract attention of the reader from the aspects 
of the physics we are going to concentrate upon, we'll con- 
sider a toy model of graphene, with free electrons being 
described by the 2d Dirac Hamiltonian. Thus the model 
can possess only a single Dirac point, and we'll present 
the existing approaches as if they were applied to this 
toy model. Notice, that due to the isotropy of the model 

x(iJ)=x(R), (2) 
where R is the distance between the sites i and j. 

APPROACHES BASED ON ZTGF 

The approach, used in Refs. [l], d, Q is based on equa- 
tion 
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X(u> = 0,q)e 



iq-R 



where 



X(w = 0,q) = 2 



d 2 k n F (£ k ) - n F (£k+ q ) 
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is the Fermi dis- 



tribution function. This approach, though looking quite 
straightforward, brings with it a problem. In a model of 
infinite Dirac cones for v;(o> = 0, q) we obtain a diverg- 
ing integral. To obtain finite values from these divergent 
integrals, as it was mentioned previously, one has to im- 
plement the complicated (and to some extent arbitrary) 
cut-off procedure 0|. 

The problem can be formulated in a different way. Be- 
ing calculated in a realistic band model, with the bands 
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of finite width, x(w = 0, q) is not a universal quantity. It 
depends not only on infrared physics, but on the proper- 
ties of electron spectrum and eigenfunctions in the whole 
Brillouin zone (even for small q). 

Another approach, formulated in Ref. [2], starts from 
a well known equation for the susceptibility 



2? f°° 

X (R) = - G 2 (R,E)dE, 

K J-oo 



(5) 



where G is the retarded green's function. Here again the 
integral diverges on both limits of integration. However 
the authors changed the contour of integration, trans- 
forming the divergent integral (J5j) into the convergent 
integral along the imaginary axis (see also Ref. [20(). 
The authors also considered RKKY interaction in gapped 
graphene, when the power law decrease of the interaction 
with the distance turns into the exponential law. Actu- 
ally, in a implicit form the authors made the transition 
from ZTGF to MGF, so our final results will be very close 
to ones obtained in Ref. [2| . 
The approach, using formula 



X(r,r') = Sn(r)/SV(r') 



(6) 



and, hence, calculating electron susceptibility on the ba- 
sis of equation 

X (R) = -- Im [G 2 (i?, E)] dE, (7) 

where Ep is the Fermi energy, was first used, in appli- 
cation to graphene to the best of our knowledge, in Ref. 
7]. An advantage of this approach is that it allows to 
easily consider the case of doped graphene, the disadvan- 
tage is that the approach, like the one presented above, 
has to deal with the divergent integral (the integral with 
respect to dE diverges at the lower limit of integration) . 
Also in this case, to obtain finite values from these di- 
vergent integrals one has to implement the complicated 
(and to some extent arbitrary) cut-off procedure. 



MGF IN FREQUENCY REPRESENTATION 



Our approach will be based on equation 21 1 
1 



x(R) 



g z (R;u))du, 



(8) 



where Q(R;iv) is the MGF in frequency momentum rep- 
resentation. Dirac equation describing electrons is 



H = v{T X k x +T y k y ), 



(9) 



where the matrix r acts in the space of two sublattices. 
From Eq. © we obtain 



1 



iuj + fi — v(r x k x + rvky)' 



(10) 



where fj, is a chemical potential. From Eq. (TTUJ) we obtain 

M 



Q CC {R^) 



kdk 



-i(u> — ifi) 
(2tt) 2 J q (w-i/i) 2 + v 2 k 2 



2w e ikR C o 3 e de = ~ ^) f°° kJ (kR)dk 
„ 2ir J (uj - ifi) 2 + v 2 k 2 

i{uj — i/i) 



2ttv 2 



-K [sign(w)(w - in)R/v] 



(11) 



Q ab {R,oj) 



k 2 dk 



(2tt) 2 J (u - i^i) 2 + v 2 k 2 



ikR cos Q-\-iQ 



de 



v f°° k 2 J l {kR)dk 



2tt J q (lj — i/i) 2 + v 2 k 2 

sign(w)(o; — ifx) 



2irv 2 



-Ki [sign(w)(w - ifi)R/v] . (12) 



[CC can mean either AA or BB; Kq and K\ are the 
modified Bessel function of zero and first order respec- 
tively). We have used mathematical identity, valid for 
Re z > [H, 



{X 2 + Z 2y j ^) dx = 2P -i r(p) K v - p+1 {cz). 



(13) 



UNDOPED GRAPHENE 



Consider first the case of undoped graphene (fj, = 0). 



MGF in frequency representation 



Using another mathematical identity [2c 



x a Kp{cx)K v {cx)dx = — 

c a l (a) 
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a + /j, + v 
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2 J \ 2 
from Eq. © we obtain [1] 



6AttvR 3 ' 
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(14) 



.(15) 



MGF in time representation 

For completeness, we reproduce here the calculation 
of x based on MGF in coordinate - imaginary time rep- 
resentation, presented in our previous paper [llj. The 
susceptibility was written as [j, [ll[ 

/•OO 

X (R) = -2 g(R;r)g(R;-T)dr. (16) 
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Transition from frequency to imaginary time representa- 
tion yields the MGF 



(17) 



As a result we obtain 

„ rr , , sign(r) '°° ' a * 
8^ 2 jo Jo 



dkk / ^ e ^cos9-^fc|r| 



G AB (R,T) = ^r dkk d e ikRc O B$+iB-vk\r\ m 

87r Jo Jo 

(18) 

Performing the angular integrations in Eq. (|18[) we get 

dkkJ (kR)e- vk \ T \ 



G cc (R;r)- SigniT) 



4/T 



1 Z" 00 

gA B (R ]T )= dkkJ 1 {kR)e~ vk ^. (19) 

4tt Jo 

(Jo and J i are the Bessel function of zero and first order 
respectively). Using mathematical identity [23| 



c n - l e- px J v {cx)dx 



(20) 



(-1) 



gn-l (vV + c 2 -p) 

lr -1 ^^ r ^ 7 



integrals in the RHS of Eq. (fTO|) can be calculated ex- 
actly, giving a well known result || 



S cc (i?;r) = 



<? as or;t) 



1 DT 

4^ (u 2 r 2 + R 2 ) 3 / 2 
1 g 
4^ (v 2 r 2 + i2 2 ) 3 / 2 ' 



(21) 



The remaining integration in Eq. (fTfj|) is trivial; as a 
result we recover Eq. (fl~5|) . 



DOPED GRAPHENE 

In the case of doped graphene the susceptibility © 
will be expressed through the integrals 

Rej^ (u-itfK$ A [(w-iv)R/v]du)^. (22) 

Considering integrals in the complex plane it is conve- 
nient to deform the contour of integration and present 
the integrals as 



'Oo />0 

Re\ I uj 2 Kl 1 {ujR/v)duj + J oj 2 Kl 1 (ujR/v)duj 



(23) 



Taking into account the identity 



K a {-ix) = -i a+1 [J a {x)+iY a {x)], 



(24) 



we get [12| 

x° c (R) = x%{R) 



xt B (R)=x A *o{R) 



k F R 

1 - III / dzz 2 J (z)Y (z) 



1G 



o 

k F R 



dzz 2 J 1 {z)Y 1 (z) 







(25) 



where kp — (m/v. The integrals in Eq. (|25[) can be pre- 
sented in terms of Meijer functions [12| (I send the reader 
to that Reference for the details). 

It is interesting to compare the RKKY exchange in 
doped graphene, with its two sublattices and linear dis- 
persion law, with that in ordinary two-dimensional elec- 
tron gas. For the latter the Green's function is 



1 



icj + /i — k 2 /2m 
Hence the susceptibility turns out to be [24 



1 



k F R. 



X(R) ~Jpj Q dzzJ Q (z)Y (z). 



(26) 



(27) 



Amusing, that the authors of Ref. [24[ were affiliated 
with the same University, as the author of the present 
paper. 



GAPPED GRAPHENE 

Consider now graphene with the gap in electron spec- 
trum described by Dirac Hamiltonian 



H = v{T X k x +T V k y )+AT Z 

The Green's function is 



(28) 



an,, a -iu-&T*-v(T*k x + Tyk y ) 

GiKuj) = u, 2 + A 2 +, 2 fc 2 ■ (29) 



MGF in frequency representation 
From Eq. (f2"9"j) we obtain 



Q UU (R,UJ): 

-iui ± A 

27TU 2 



-ioj =f A f°° kJ {kR)dk 
2tt J q uj 2 + A 2 +v 2 k 2 



K (V^ 2 + A 2 i?/v) 



(30) 



AB -v r ej^dk 

y {RiU) = — 



-V^ 2 + A 



2ttv 



2tt J u 2 + A 2 +v 2 k 2 
, K x (V^ 2 + A 2 i?/ V ) . (31) 
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In Eq. minus corresponds to AA and plus to BB. 

Substituting into Eq. ([5]) we obtain [25| 



x cc {R) 



47T 3 ?; 4 

{uj 2 - A 2 )K 2 (v^+tfiJ/v 



dui. 



(32) 



x AB {R) = 



l 



47T 3 U 4 

(w 2 + A 2 )^ 2 (V ^ 2 + A 2 i?/f) dw. (33) 



The remaining integrations can be performed analyti- 
cally in two limiting cases. 

Consider first the case AR/v -C 1. Here it is appro- 
priate to mention the relation between the toy model, 
we are using, and real graphcnc. The existence of two 
Dirac points in graphene leads to additional angular de- 
pendent factor in the formula for the RKKY interaction. 
It was thoroughly studied previously and does not in- 
terfere with the physics we are discussing in this work. 
More interesting is the condition of the applicability of 
the infinite Dirac cones dispersion law to the calculation 
of RKKY interaction in graphene. From Eq. (|2~lj) we 
see that characteristic values of r entering into integral 
(HHJ) are of the order of R/v; hence characteristic k en- 
tering into integrals (TiH|) are of the order of 1/R. Thus 
the approximation of linear dispersion law is applicable, 
provided R ^> a, where a is the graphene lattice con- 
stant. Now we realize, that the case AR/v <C 1 can be 
described in the framework of the model if A <C av, or 
in simple terms, if the gap is narrow in comparison with 
the graphene band width, which is certainly true in most 
cases of gapped graphene. So finally, in the case of a 
narrow gap and relatively small distances we can go to 
the limit A — > in Eqs. (|32p and Q33p. and recover the 
results of the gapless case. 

The other limiting case AR/v 3> 1 is more interest- 
ing. In this case we may use asymptotic expression for 
modified Bessel functions 



K v {z) 



(34) 



After calculating the resulting integrals in Eq. (f3"3"|) using 
the Laplace method, we obtain the same result both for 
inter-sublattice and intra-sublattice susceptibility 



1 



X 



8v V ttvR 



3/2 



-2RA/v 



(35) 



It is worth paying attention to the fact that Eq. 
seems to contradict rigorously proved theorem stating 
that for any bipartite lattice at half filling, the RKKY in- 
teraction is antifcrromagnetic between impurities sitting 
on top of atoms belonging to opposite sublattices (i.e., 
A and B sublattices in graphene), and is ferromagnetic 



between impurities sitting on top of atoms belonging to 
the same sublattice [1, 11, (26|. However, the theorem 
is not applicable to Hamiltonian (|28|) , with its last term 
meaning that if we rewrite the Hamiltonian in the tight- 
binding representation, the intra-sublattice hopping will 
appear, hence the lattice is no longer bipartite. More 
specifically, the spectrum still has the symmetry of that 
in bipartite lattice, but the wave functions do not 11|. 



MGF in time representation 

Here we'll restrict ourselves with the calculation of 
X AB ■ Transition from frequency to imaginary time rep- 
resentation yields the MGF 

1 vk 



G AB (k,r) = ~ 



2 VA 2 + v 2 k 2 
Hence instead of Eq. (TiU]) we obtain 

J (kR)e-^+^^ 

_P 

VA 2 +v 2 k 2 ' 

Using mathematical identity [2S 



(36) 



(37) 



g^(fl;r) = f / dk *" J^e-^^H. 
4tt J VA 2 + v 2 k 2 

(38) 



oV x 2 +z 2 

\Jx 2 + z'- 



-J v (cx)dx = B v , (39) 



where 



B x = c 1 + zy/p 2 + ■ 



(J5 2 + C 2 ) 3 / 2 



we obtain 



^-Vv 2 t 2 +R 2 A/v 



B» W r) . |(l + V^TWNv) + 



(40) 



Thus we obtain 



,AB 



732 POO ■ 

(R)=-—J (l + Vv 2 T 2 +R 2 A/vj 



-2Vv 2 t 2 +R 2 A/v 
(V 2 T 2 + R 2 f 



-dr. 



(41) 



In the case AR/v > 1 we can calculate integral in Eq. 
(|4"Tj) using the Laplace method, to recover Eq. ([33)1 . 



CONCLUSIONS 

In the end we would like to mention again that though 
we were considering the case of T = 0, we have found, 
that, as it is not infrequently happens, the MGF have 
advantages over ZTGF. In particular, using the former 
one have to operate only with the convergent integrals, 
in distinction to what happens when one uses the latter. 
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